Half-metallicity and magnetism of GeTe doped with transition metals V, Cr and Mn: 
a theoretical study from the viewpoint of application in spintronics 
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This work presents results for the magnetic properties of the compound GeTe doped with 3d 
transition metals V, Cr and Mn from the viewpoint of potential application in spintronics. We report 
a systematic density-functional study of the electronic structure, magnetic and cohesive properties 
of these ternary compounds in both rock salt and zinc blende structures. In both cases, it is the 
Ge sublattice that is doped with the three transition metals. Some of these compounds are found 
to be half-metallic at their optimized cell volumes. For these particular cases, we calculate both 
exchange interactions and the Curie temperatures in order to provide some theoretical guidance to 
experimentalists trying to fabricate materials suitable for spintronic devices. Discussions relating 
our results to the existing experimental studies are provided whenever applicable and appropriate. 
Apparent discrepancy between experimental observations and our theoretical result for the case of 
Mn-doping is discussed in detail, pointing out various physical reasons and possible resolutions of 
the apparent discrepancy. 

PACS numbers: 71.20.Nr, 71.20.Lp, 75.30.Et, 75.10.Hk 
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I. INTRODUCTION 

In recent years GeTe thin films doped with 3d tran- 
sition metals have received considerable attention from 
experimentalists interested in the field of magnetic 
semiconductors >*— Ferromagnetic order was observed for 
the Cr, Mn, and Fe doped films, whereas Ti, V, Co and Ni 
doped films were found to be paramagnetic^ The Curie 
temperatures T c of these thin films have been found to 
depend on the type and concentration of the dopants, 
with a high value of 140 K, reported for Gei-^Mn^Te 
for x=0.51»i Later on this group of researchers reported 
even higher value of T c ~ 190 K in dilute magnetic semi- 
conductor (DMS) Ge . 92 Mn . 08 Te.. 7 T c s around 200-250 
K have been reported by Fukuma et al£ for Gei-^Cr^Te 
thin films for low values of x (< 0.1). The thin films in 
all these studies were of predominantly rock salt struc- 
ture. For spintronic applications, the desired proper- 
ties of such ferromagnets are half-metallicity and rel- 
atively high Curie temperatures. Possibility of half- 
metallicity in Cr- and V-substituted GeTe bulk com- 
pounds has been reported by Zhao et air- based on the- 
oretical calculations using linear augmented plane waves 
(LAPW) WIEN2k code for some ordered Ge-V-Te and 
Ge-Cr-Te compounds. These authors left the issue of 
the Curie temperature unattended, as only the electronic 
band structure, density of states and charge density were 
studied, and not the exchange interaction in these com- 
pounds. A theoretical study of the electronic structure 
of Gei-zMn-rTe, based on the generalized gradient ap- 



proximation plus Hubbard U (GGA+U), has been pre- 
sented by Ciucivara et al. 10 Among the magnetic proper- 
ties only the magnetic moment was considered, the issue 
of the exchange interaction and T c was not addressed. 
Recently a large number of experimental studies related 
to magnetism and magnetic transitions in several Ge-Te 
based alloys, such as Ge-Cr-Te, Ge-Mn-Te and Ge-Fe-Te) 
have been reportedf^— In addition, there is substantial 
amount work on DMSs in general (for a review see [22|h 
where the magnetic properties arise via percolation ef- 
fects among a very small concentration of the magnetic 
atoms. For a discussion on magnetic percolation in DMSs 
readers may consult the work by Bergqvist et al. 23 

In view of the large body of existing experimental work 
as mentioned above, we have embarked on a systematic 
theoretical study of the GeTe system doped with 3d tran- 
sition metal atoms. The most commonly occurring struc- 
ture for both bulk and thin films of these compounds is 
the NaCl or the rock-salt (RS) structure. However, from 
the theoretical viewpoint an equally interesting structure 
to study is the zinc blende (ZB) structure. Both RS and 
ZB are fec-based, but differ in terms of the distance be- 
tween the Ge and Te atoms. Thus irrespective of whether 
one dopes the Te- or the Ge-sublattice with a magnetic 
(3d) atom, the magnetic effects are expected to differ be- 
cause of different levels of hybridization between the 3d 
orbitals and the s and p orbitals of the Ge or Te atoms. 
Experimentally it may be possible to grow both RS and 
ZB structures, even though the ground state structure 
appears to be RS. Thus we have studied the electronic 
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and magnetic properties of both these structures for dif- 
ferent concentrations of the magnetic V, Cr, and Mn 
atoms. We have employed supercell method as well as 
the coherent potential approximation (CPA) to study the 
effect of doping at various concentrations. The supercell 
calculations are carried out using the full potential lin- 
ear augmented plane wave (FP-LAPW) method using 
the WIEN2k code^i The CPA calculations are carried 
out within the frame-work of the linear muffin-tin orbital 
(LMTO) method using the atomic sphere approximation 
(ASA) ; 25 ' 26 In a limited number of cases, when the half- 
metallicity is predicted on the basis of both supercell and 
CPA calculations, we have provided results for exchange 
interactions and Curie temperatures for these alloys. 



II. COMPUTATIONAL DETAILS 

The crystal structures of the ternary GeTe based com- 
pounds Gei-zTM^Te (TM stands for transition metal 
atoms V, Cr and Mn) were constructed from the unit 
cell of RS and ZB structures as follows. The doping lev- 
els of x=0.25 or 0.75 were achieved by replacing the Ge 
atoms at the vertex site or face-center sites, respectively, 
of the RS/ZB unit cell , with the TM atoms. Both cases 
(x=0.25 or 0.75) have the same space group (Pm3m, or 
221) for the RS structure, and (P43m, or 215) for ZB. 
The x=0.5 doping is obtained by replacing the Ge atoms 
at the four compatible face-center sites. In this way we 
generate tetragonal structure (P4m2, 115) for the ZB 
case and (P4/mmm, 123) for RS. This procedure of real- 
izing the doping levels results in minimal primitive cells, 
with the highest possible symmetry. Numerous semicon- 
ductors are known to crystallize in the ZB or RS struc- 
tures. As such these ternary compounds should be com- 
patible with a large number of semiconductors. 

Calculations were carried out within the framework 
of the density- functional theory (DFT),— using the 
WIEN2k 24 code based on the full-potential linear aug- 
mented plane wave plus local orbitals method. The 
generalized gradient approximation (GGA) proposed by 
Perdew et al. was used for exchange and correlation 
potentials. 28 ' 29 We consider full rclativistic effects for the 
core states and use the scalar-relativistic approximation 
for the valence states, neglecting the spin-orbit coupling. 
The latter is known to produce only a small effect on the 
results that are of interest in this work, e.g., density of 
states and the energy gaps. 

As a note to practitioners of this code, we used 3000 
K points (Monkhorst-Park grid}2£ for the Brillouin-zone 
integration, set Rmt * Kmax = 8.0 and carried out the 
angular momentum expansion up to l ma x = 10 in the 
muffin tins, and used G max — 14 for the charge density. 
All core and valence states are treated self-consistcntly. 
Convergence with respect to basis set and k-point sam- 
pling was carefully verified. The self-consistent calcula- 
tion is allowed to stop only when the integrated charge 
difference per formula unit, J \p n — p n -i\dr, for the in- 



put and output charge densities p n -\ and p n is less than 
10~ 4 . In the calculation of Ge-X-Tc in different struc- 
tures, the muffin-tin (MT) radii are chosen to be 2.3, 2.2 
and 2.5 a.u. for Ge, X (X=V, Cr, Mn) and Te atoms, 
respectively. 



III. ELECTRONIC STRUCTURE 

A. Results of supercell calculations using the 
FP-LAPW method 

For each case equilibrium lattice parameter was ob- 
tained by minimizing the total energy with respect to 
the cell volume. All electronic properties such as the 
density of states (DOS), energy bands and magnetic mo- 
ments were then calculated for the equilibrium lattice 
parameters. Among all the ternary TM compounds with 
the doping levels considered, we find 9 half-metallic (HM) 
cases: one Mn-doped, three Cr-doped, and three V-doped 
cases for the RS structure; and one Cr-doped and one 
V-doped cases for the ZB structure (see Tables I-IV). 
There is a small drop in the value of the equilibrium lat- 
tice parameter with increasing dopant concentration in 
all cases. This is supported by the measurements of lat- 
tice parameters for the thin films of Gei^Mn^Te and 
Gei-zCr^Te grown in RS structure^ 5 - Bulk moduli, cal- 
culated by using Birch-Murnaghan equation of stated 
are found to increase with the doping level for V- and 
Cr-doping. For Mn-doping this change is non-monotonic. 
The equilibrium lattice constants, X-Te bond lengths, 
and bulk moduli, magnetic moments, minority-spin gaps, 
and half-metallic gaps for the ternary compounds in RS 
structure are summarized in Tables I and II. The same 
quantities for the ternary compounds in ZB structure are 
summarized in Tables III and IV. In general, the band 
gap tends to increase with decreasing lattice parameter. 
This is understandable as the increased hybridization due 
to decreasing inter-atomic distances leads to increased 
separation between the energies of the bonding and anti- 
bonding states. 

Tables [II] and IIVI show that most of the moments are 
associated with the muffin-tins around the TM atoms 
and the interstitial space. The orbitals associated with 
Ge and Te atoms have little spin polarization. While the 
values of the TM and interstitial moments are dependent 
on the muffin-tin radii, the integer values of the total 
moments shown in the last column clearly identify the 
half-metallic cases. 

Note that in all cases studied and reported in Tables 
Hl lIVI the RS structure has lower energy than the ZB, 
indicating the equilibrium bulk phase to be RS at low 
temperatures. However, it might be possible to grow thin 
films of these compounds in ZB structure under suitable 
conditions. Therefore, there is some merit in comparing 
the magnetic properties of these two fec-based phases. 
So far all thin films of these compounds seem to have 
been grown in RS structure or with small rhombohedral 
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TABLE I: Results obtained via the WIEN2k code for 
the RS compounds Ge 4 _„X„Te 4 (X=V,Cr,Mn) : the 
equilibrium lattice constants a, the X-Te bond length 
LxTe, bulk modulus B, minority-spin gaps Gmis (eV), 
and half-metallic gaps Ghm (eV). All results shown are 
obtained using GGA (see text). 



compounds 


o(A) 


L(A) 


B(GPa) 


G M is(eV) 


G HM (eV) 


Ge 3 ViTe 4 


5.9283 


2.9642 


51.3 


0.599 


0.097 


Ge 2 V 2 Te 4 


5.8497 


2.9249 


55.8 


0.816 


0.170 


GeiV 3 Te 4 


5.7865 


2.8933 


59.7 


1.007 


0.213 


Ge 3 CriTe 4 


5.9325 


2.9662 


51.4 


0.571 


0.220 


Ge 2 Cr 2 Te 4 


5.8579 


2.9290 


54.9 


0.789 


0.382 


GeiCr 3 Te 4 


5.7860 


2.8930 


58.2 


0.980 


0.371 


Ge 3 MmTe 4 


5.9638 


2.9819 


49.7 


0.710 


0.200 


Ge 2 Mn 2 Te 4 


5.9139 


2.9568 


46.1 






GeiMn 3 Te 4 


5.8413 


2.9206 


43.1 







TABLE III: Results obtained via the WIEN2k code for 
the ZB compounds Ge 4 _„X„Te 4 (X=V,Cr,Mn) : the 
equilibrium lattice constants a, the X-Te bond length 
LxTe, bulk modulus B, minority-spin gaps Gmis (eV), 
and half-metallic gaps Ghm (eV). All results shown are 
obtained using GGA (see text). 



compounds 


a(A) 


L(A) 


B(GPa) 


Gmis (cV) 


G HM {eV) 


Ge 3 ViTe 4 


6.5842 


2.8510 


27.7 






Ge 2 V 2 Te 4 


6.4316 


2.7953 


33.8 






GeiV 3 Te 4 


6.3109 


2.7327 


42.2 


1.061 


0.482 


Ge 3 CriTe 4 


6.5950 


2.8556 


30.4 






Ge 2 Cr 2 Te 4 


6.4347 


2.7866 


34.0 






GeiCr 3 Te 4 


6.3280 


2.7401 


38.8 


1.144 


0.253 


Ge 3 MmTe 4 


6.6405 


2.8754 


29.1 






Ge 2 Mn 2 Te 4 


6.5305 


2.8281 


30.1 






GeiMn 3 Te 4 


6.4222 


2.7809 


26.5 







TABLE II: Magnetic moment per magnetic atom (f.u.) 
in units of bohr magneton /is, for the RS compounds 

Ge 4 _„X„Te 4 (X=V, Cr, Mn), at their respective 
equilibrium volumes, obtained by using the FP-LAPW 
method (WIEN2k). Moments associated with X, Te, Ge 
muffin-tin spheres and the interstitial region (Int) are 
shown separately. 



compounds X(/is) Te(/is) Ge(/zs) Int(/is) Tot (/is) 



Ge 3 ViTe 4 


2.463 


-0.050 


0.026 


0.624 


3.000 


Ge 2 V 2 Te 4 


2.482 


-0.053 


0.043 


0.606 


3.000 


GeiV 3 Te 4 


2.497 


-0.076 


0.048 


0.589 


3.000 


Ge 3 CriTe 4 


3.580 


-0.065 


0.017 


0.576 


4.000 


Ge 2 Cr 2 Te 4 


3.575 


-0.065 


0.031 


0.565 


4.000 


GeiCr 3 Te 4 


3.567 


-0.120 


0.047 


0.558 


4.000 


Ge 3 MmTe 4 


4.243 


0.035 


0.024 


0.583 


5.000 


Ge 2 Mn 2 Te 4 


4.226 


0.032 


0.053 


0.561 


4.947 


GeiMn 3 Te 4 


4.140 


0.050 


0.083 


0.494 


4.167 



TABLE IV: Magnetic moment per magnetic atom 

(f.u.) in units of bohr magneton [ib, for the ZB 
compounds Ge 4 _„X„Te 4 (X=V, Cr, Mn), at their 
respective equilibrium volumes, obtained by using the 
FP-LAPW method (WIEN2k). Moments associated 
with X, Te, Ge muffin-tin spheres and the interstitial 
region (Int) are shown separately. 



compounds X (/is) Te(/is) Ge(/is) Int(/is) Tot(/is) 



Ge 3 ViTe 4 


2.535 


-0.079 


0.022 


0.719 


3.007 


Ge 2 V 2 Te 4 


2.510 


-0.109 


0.027 


0.690 


3.009 


GeiV 3 Te 4 


2.463 


-0.113 


0.040 


0.674 


3.000 


Ge 3 CriTe 4 


3.706 


-0.079 


0.013 


0.675 


4.103 


Ge 2 Cr 2 Te 4 


3.625 


-0.126 


0.014 


0.658 


4.045 


GeiCr 3 Te 4 


3.560 


-0.145 


0.021 


0.627 


4.000 


Ge 3 MmTe 4 


4.171 


0.024 


0.033 


0.634 


5.003 


Ge 2 Mn 2 Te 4 


4.190 


0.054 


0.068 


0.635 


5.002 


GeiMn 3 Te 4 


4.191 


0.056 


0.078 


0.597 


4.889 



distortions from this structure ,ii2iLiii 

The densities of states for some these ordered com- 
pounds at their equilibrium lattice parameters are shown 
in Figs.([T][2]). These results clearly show in which cases 
half-mctallicity is most robust, i.e. the Fermi level is 
most widely separated from the band edges. The partial 
atom-projected densities of states shown for Ge2V2Te 4 , 
shows that the states at the Fermi level have a high TM 
character. The situation should be similar for TM=Cr, 
Mn as well, with the amount of the TM character chang- 
ing with TM concentration. The 25% Mn-doping case 
needs special attention, since the metallic character is ei- 
ther marginal (ZB structure) or absent (RS structure). 
For Gc-Tc in the ZB structure, 25% Mn-doping of the 
Ge-sublattice produces a borderline, i.e., zero gap, semi- 
conductor or semimetal. For the RS-structure, the same 



Mn-doping level of the Ge-sublattice produces a narrow 
gap semiconductor. As we will see in the next section, 
the nature of magnetism in these compounds is dictated 
by the availability of free carriers or holes. Thus the Mn- 
doping case will be considered in some additional detail. 

The results presented in Tables I-IV and Figs. 1-2 
are for ordered alloys, while the TM-doped thin films 
of GeTe are partially random, in the sense that the 
Ge-sublattice sites are occupied randomly by Ge and 
TM atoms. We have studied the electronic struc- 
ture of these partially random alloys using the tight- 
binding linear muffin-tin orbital(TB-LMTO) method in 
conjunction with the CPA i 25 ' 26 These calculations em- 
ployed exchange-correlation potential of Vosko, Wilk and 
Nusair, 32 an s,p,d,f basis set, relativistic treatment of 
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FIG. 1: Densities of states in Ge 3 TMiTe 4 (TM=V,Cr,Mn) alloys in RS and ZB structures. 
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FIG. 2: (Color online) Total and component-resolved 
densities of states in V2Ge2Te4 alloys in RS and ZB 
structures. 



core electrons and scalar-relativistic treatment of valence 
electrons. The results, particularly with respect to half- 
metallicity, were similar to those of the ordered alloys, 
apart from expected smoothing of some peaks in the 
DOS. The gaps values were marginally lower and can be 



ascribed to the differences between local density approx- 
imation (LDA) in TB-LMTO and GGA in FP-LAPW. 
The equilibrium lattice parameters were 2-4% higher in 
TB-LMTO LDA calculations. 



IV. EXCHANGE INTERACTION AND CURIE 
TEMPERATURE 

First-principles calculations of the thermody- 
namic properties of itinerant magnetic systems, via 
mappin g 33 ' 34 of the zero temperature band energy onto 
a classical Heisenberg model: 

H c g - 



1^ Ji i ' 

i,3 



(1) 



has been discussed in detail in our previous 
publications! 35 ' 36 Here i,j are site indices, e, is the 
unit vector pointing along the direction of the local 
magnetic moment at site i, and Jjj is the exchange 
interaction between the moments at sites i and j. The 
calculations are based on a mapping procedure due to 
Liechtenstein et ali^r— The method was later extended 
to random magnetic systems by Turek et al. , using 
CPA and the TB-LMTO method^ The exchange 
integral in Eq.([T]) is given by 



Jn = — lim Im I tTj, 

47T e->0+ 



AiWgUzfaWgi dz, (2) 



where z = E + ie represents the complex energy variable, 
L = (l,m), and Aj(z) = P}{z) — Pf{z), representing 
the difference in the potential functions for the up and 
down spin electrons at site V. gfj(z)(a =t, I) repre- 
sents the matrix elements of the Green's function of the 
medium for the up and down spin electrons. For sublat- 
tices with disorder, this is a configurationally averaged 
Green's function, obtained via using the prescription of 
CPA. It should be noted that the spin magnetic moments 
are included in the above definition of Jij. Positive and 
negative values of Jy imply FM and AFM couplings, re- 
spectively, between the atoms at sites i and j. 

Existence of the HM gap in spin polarized electronic 
densities of states does not guarantee that the substance 
is actually ferromagnetic (FM). The search for the ground 
magnetic state often needs to be guided by ab initio 
calculation of the exchange interactions between various 
atoms. In the above procedure outlined by Liechtenstein 
et al.^r— exchange interactions are calculated by con- 
sidering spin deviation from a reference state. Negative 
exchange interactions resulting from calculations based 
on a FM reference state would suggest instability of the 
assumed FM ground state. With this mind, we have com- 
puted the exchange interactions for the most interesting 
cases, i.e. those promising robust HM states. The results 
for the exchange interactions are shown in Figs.([5][5]) 

In Figs. [3]and[Hwe show results for the cases with 75% 
doping with Cr and V for the RS and ZB structures, as 
this level of doping yields largest values of the minor- 
ity spin and HM gaps. In Fig. [5] we consider the 25% 
Mn-doping case for the RS structure, the only Mn-doped 
case studied which shows the promise of half-metallicity. 
The solid lines (with circles) in these figures refer to the 
alloys with the equilibrium lattice parameters. In or- 
der to understand some trends, we have also considered 
lattice parameters above and/or below the equilibrium 
values. Ferromagnetic interactions are strongest in the 
Cr-doped ZB GeTe. The Cr-doping case is most promis- 
ing, as in both RS and ZB structures FM interactions 
dominate, and these remain FM with changes in the lat- 
tice parameter. The next promising case is V-doping in 
the ZB structure, the RS counterpart showing strong an- 
tiferromagnetic (AFM) interactions, particularly at and 
around the equilibrium lattice parameter. These AFM 
interactions weaken on both sides of the equilibrium lat- 
tice parameter value and become FM only at much higher 
volume. FigJS] shows that all interactions up to many 
neighbor shells are AFM, and then die off to zero. Hence 
the ground state cannot be FM. This conclusion needs to 
reinforced by examining the Lattice Fourier transforma- 
tion of all the interactions. However, the preponderance 
of AFM interactions would dictate that the ground state 
is perhaps AFM or of complex magnetic structure. This 
issue is further explored in the next section. For the 25% 
Mn-doped solid the interactions are strongly AFM not 
only at the equilibrium lattice parameter, but also for 
a substantial range of the lattice parameter around the 
equilibrium value. The interactions can become FM only 
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FIG. 3: (Color online) Exchange interaction between 
the V atoms in Geo.25Vo.75Te as a function of 
interatomic distance expressed in units of lattice 
parameter a. Solid lines with solid circles show the 
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FIG. 4: (Color online) Exchange interaction between 
the Cr atoms in Geo.25Cro.75Te as a function of 
interatomic distance expressed in units of lattice 
parameter a. Solid lines with solid circles show the 
results for the equilibrium lattice parameters. 



at unrealistically large lattice parameter, while compres- 
sion would result in stronger AFM interactions. While 
our theoretical calculations suggest that 25% Mn-doping 
of the Ge-sublattice should lead to a ground state that 
is either AFM or of a more complex magnetic structure, 
some recent experimental results suggest ferromagnetism 
in these or similar Mn-doped compounds. Thus we will 
discuss the Mn-doping case separately in order to shed 
some light on the apparent discrepancy between theo- 
retical results and the available experimental studies. In 
the following paragraphs we first present estimates of the 
Curie temperature for the Cr- and V-doped Ge-Te. 

We have calculated the Curie temperature T c using 
both the mean-field approximation (MFA) and the more 
accurate random-phase approximation (RPA). 41 If the 
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FIG. 5: (Color online) Exchange interaction between 
the Mn atoms in Geo.75Mno.25Te as a function of 
interatomic distance expressed in units of lattice 
parameter a. Solid lines with solid circles show the 
results for the equilibrium lattice parameters. 
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FIG. 6: (Color online) T c (via RPA together with 
VCA) versus lattice parameter a in Gei-^X^Te (X=V, 

Cr) alloys in ZB structure. The arrows indicate the 
values for the equilibrium lattice parameters for £=0.75 



magnetic sublattice consists only of the magnetic atoms 
X, then in the MFA, the Curie temperature is given by 



, rpMFA 2 \ - x,X / Q \ 
i#0 

where the sum extends over all the neighboring shells 
and involves the exchange interactions between the mag- 
netic atoms X. MFA is known to grossly overestimate T c . 
A much more improved description of finite-temperature 
magnetism is actually provided by the RPA. Again, if the 
magnetic sublattice consists only of the magnetic atoms 
X, then the RPA T c given by 

EtJ^w-j^q)]- 1 . (4) 

q 

Here N denotes the order of the translational group ap- 
plied and J x,x (q) is the lattice Fourier transform of the 
real-space exchange interactions </*' X - In order to ad- 
dress the randomness in the Ge-TM sublattice, we have 
modified Eqs.Q and (|4]) for our Gci-^X^Te alloys using 
virtual crystal approximation (VCA). This involves sim- 
ply weighting the exchange integrals in Eq. (fTJ) by x 2 , 
where x is the concentration of the X (TM) atoms. As a 
result, Tc's obtained from Eqs.([5]) and ^ get multiplied 
by the same factor. In this way the problem is formally 
reduced to a nonrandom case. This approximation fails 
for low concentrations, below the percolation limits The 
error decreases monotonically for higher and higher con- 
centrations. In general, the VCA results may somewhat 
overestimate the Curie temperature. 

A problem arises in the computation of T c using cither 
Eq.Q or (jlj when, in addition to the robust moments 
on the magnetic atoms, there are induced moments on 



apparently non-magnetic atoms, interstitial spaces or, in 
case of the LMTO method, empty spheres. This problem 
has been discussed in detail in our previous publication 35 
and references cited therein. As shown by Sandratskii et 
alM- the calculation of T c using RPA is considerably more 
involved even for the case where only one secondary in- 
duced interaction needs to be considered, in addition to 
the principal interaction between the strong moments. 
The complexity of the problem increases even for MFA, 
if more than one secondary interaction is to be consid- 
ered. In our case, the induced moments are small and as 
such interactions involving non- magnetic atoms/spheres 
can be neglected in the first approximation. In addition, 
we are only interested in a rough estimation of T c , the ob- 
ject being able to determine which doping would lead to 
a higher T c and hopefully be close to or above the room 
temperature. Our results, obtained by using the RPA 
and considering only the TM moments in Eq. @ , and fur- 
ther modified in the spirit of VCA, are summarized Figs. 
Oand \7\ We show T c as a function of the lattice parame- 
ter for various relevant levels of doping. The equilibrium 
lattice parameter values are indicated with arrows. Dop- 
ing with Cr or V in the ZB structure for doping levels 
around 75% is promising, with V-doping showing higher 
T c values than Cr-doping. V-doping in RS structure is 
not recommended, while Cr-doping in RS structure leads 
to reasonably high values of T c for all levels of doping 
higher than 50%. In a previous publication 3 ^ we have 
already presented a comparison of pure CrTe in ZB and 
RS structures, pointing out various advantages of con- 
sidering RS CrTe over ZB CrTe from the viewpoint of 
spintronics applications. 
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FIG. 7: (Color online) T c (via RPA together with 
VCA) versus lattice parameter a in Gei-^X^Te (X=V, 

Cr) alloys in RS structure. The arrows indicate the 
values for the equilibrium lattice parameters for x=0.75, 
and 0.50 for Ge^Cr^Te and x=0.75 for Gei_ x V x Te. 



V. THE CASE OF Ge-Mn-Te 

Our results for the alloys Ge3Mn!Te 4 and 
Geo.75Mno.25Te need to be examined carefully in 
the context of the existing experimental studies, many 
of which claim evidence of ferromagnetism over a wide 
range of Mn concentration with T c 's up to 100-200 K. 
Our LMTO-ASA-based exchange interactions indicate 
that the ground state should not be FM, at least for the 
25% Mn concentration case studied by us. These results 
are based on LDA. So first we have checked that the use 
of GGA and more accurate full-potential results are not 
drastically different. In Table [V] we show the results for 
the FP-LAPW (Wien2k) results obtained with GGA for 
the ordered Ge3MniTe4 alloy. These results show that 
for both RS and ZB structures the AFM state energies 
are lower than the FM state energies, with the AFM 
[111] configuration having the lowest energy among the 
three AFM configurations considered. The AFM [111] 
unit cells of Ge3MniTe4, doubled with respect to the 
conventional FM unit cells, are shown in Fig. [8] 

In the RS structure the energy difference of the AFM 
configurations, in particular the [111] case, is larger than 
for the ZB case. To examine the magnetic structure of the 
ground state we also consider the lattice Fourier trans- 
form of the exchange interaction between the Mn atoms: 

J (°0 = H j or :Mn ex P (*q • R ) ; ( 5 ) 

R 

where q is a wave vector in the BZ of the fee lattice. A 
maximum at the L point would point to the ground state 
being AFM [111], while that at the T point would indi- 
cate the ground state being FM. We have examined the 
ground state magnetic structure via J(q) using both FM 
and the disordered local moment (DLM)i2^2 reference 
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TABLE V: Energy differences between ferromagnetic 

(FM) and various antiferromagnetic (AFM) 
configurations of RS and ZB GeaMniTe^ The energies 
of the AFM unit cells are compared with 2 FM unit 
cells. MO stands for magnetic order. The energy 
difference AE of the AFM configurations with respect to 
the FM, in the last column, are given per formula unit. 



MO 


Energy (eV) 


AE (meV) 


RS structure: 






FM 


-942305.912 





AFM001 


-1884611.922 


-49 


AFM110 


-1884611.981 


-79 


AFM111 


-1884612.033 


-105 


ZB Structure: 






FM 


-942304.475 





AFM001 


-1884608.979 


-14 


AFM110 


-1884609.006 


-28 


AFM111 


-1884609.018 


-33 



states. Within the Stoner model, a nonmagnetic state 
above the Curie temperature T c is characterized by the 
vanishing of the local moments in magnitude. This ob- 
viously flawed description of the nonmagnetic state can 
be improved by using the DLM model, where the lo- 
cal moments remain nonzero in magnitude above T c , but 
disorder in their direction above T c causes global mag- 
netic moment to vanish. For Mn this is known to be 
a good approximation. Within the collinear magnetic 
model, where all local axes of spin-quantization point in 
the same direction, DLM can be treated as a binary al- 
loy problem^ The results for the RS structure, obtained 
by using the FM and DLM reference states are shown in 
sections (a) and (b) of Fig. [§] respectively. Exchange in- 
teractions based on both FM and DLM reference states 
indicate that the ground state magnetic structure of RS 
Geo.75Mno.25Te is AFM [111]. 

However, the situation for ZB Geo.75Mno.25Te is not 
so clear. Both DLM and FM reference states yield J(q) 
curves with a broad maximum enclosing the symmetry 
points X, W, and K. In addition, the values of J(q) 
at these points is marginally higher than that at the L 
point (sections (a) and (b) of Fig. ITD|) . What is cer- 
tain is that the ground state is not FM, with J(q) be- 
ing a minimum at the T-point. Thus the ground state 
for ZB Geo.75Mno.25Te may involve a complex magnetic 
structure or the substance may enter a spin-glass state at 
low temperatures. In order to ensure that our results do 
not suffer from a convergence problem, we have looked 
at the J(q) values as a function of increasing shells of 
neighbors. The results for ZB Geo.75Mno.25Te with FM 
reference state are shown in section (c) of FigJTO] Results 
for 10 shells of neighbor are indistinguishable from those 
for 63 shells considered in all our calculations, indicating 
that our results are well-converged. 
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(a) RS AFM[111] (b) ZB AFM[111] 

FIG. 8: (Color online) Unit cells of Ge 3 M ni Te 4 in AFM fill] configuration. 




FIG. 9: (Color online) Lattice Fourier transform of 
Mn-Mn interactions for various lattice parameters in RS 
Geo.75Mno.25Te for (a) FM and (b) DLM reference 
states. Solid lines represent the results for the 
equilibrium lattice parameter. 



Experimentally, Gei^Mn^Te alloys, with X=Cr, 
Mn, Fe have been reported to exhibit ferromagnetic 
behavior— ~— over wide range of the concentration x. 
The structure of these compounds, grown as thin films, 
is usually found to be RS or small deviations from this. 
For Gei_j;Mn x Te the structure is found to be rhom- 
bohedrally distorted NaCl for small Mn concentrations, 
with the distortion vanishing for x > 0.18ii These stud- 
ies also report presence of carriers in the samples, with 
transport measurements most often indicating p-type 
carriers jliSilii independent of temperature. Fukuma and 
co-workers identify the carriers as holeSfii^ Depending 
on carrier concentration, brought about by varying Mn 
concentration or hydrostatic pressure, T c may vary non- 
monotonicallyJ*!^ Highest reported T c is around 130-140 
K for Mn concentration x slightly over 50% in the Ge- 




FIG. 10: (Color online) Lattice Fourier transform of 
Mn-Mn interactions in ZB Geo.75Mno.25Te for (a) FM 
and (b) DLM reference states and for various lattice 
parameters. Section (c) shows variation with respect to 
the number of neighbor shells. Results for 10 shells of 
neighbors are indistinguishable from those of 63 shells 
used in all our calculations. Solid lines in (a) and (c) 
represent the results for the equilibrium lattice 
parameter. 



Mn sublattice, and for x = 0.25 T c is ~ 70 KX Other 
authors report T c values more or less in the same range 
depending on Mn concentrationj2ii£ T c values reported 
for RS or nearly RS Gei_ x Cr K Te films are higher^ Note 
that the lattice parameters reported for x = 0.24 and 
0.55 in Gei_ x Mn x Te are 5.939 and 5.895 A, respectively, 
which compare well with our calculated values 5.964 and 
5.914 A, for x = 0.25 and 0.50, respectively (Table H]). 

Our zero temperature results indicate that exactly at 
the concentration #=0.25, Gei-zMn^Te should be either 
a semiconductor with small gap or at best a zero gap 
semiconductor, in both RS and ZB structures. The mag- 



netic structure of the ground state in this case is not 
FM, it is either AFM or more complex. At finite tem- 
peratures, free carriers in the conduction band and holes 
in the valence band are expected and the number of such 
carriers may be significant in view of the fact that the 
gap is close to zero. However, it is believed that there 
are also a large number of carriers, mainly holes, with 
temperature-independent concentration. 2 RS GeTe is re- 
ported to be a narrow gap semiconductor with a band 
gap of about 0.2 eV4^ Its very high p-type conduction 
is ascribed to large cation (Ge) vacancies ) 15 ! 44 " — which 
seem to dominate over transport due to thermally ex- 
cited carriers. We have thus decided to explore this is- 
sue, in the framework of our zero temperature formal- 
ism, in several ways. Guided by the fact that experi- 
ments indicate presence of uncompensated carriers, with 
temperature-independent concentrations, we first explore 
the effects of such carriers using the simplest possible ap- 
proach. The easiest thing for us is to find the change 
in Mn-Mn exchange interaction by moving the Fermi 
energy up (simulating holes) or down (simulating elec- 
trons). This is a one step non-self-consistent (frozen po- 
tential) calculation, in the spirit of the rigid band model. 
After self-consistency has been achieved for the equilib- 
rium lattice parameter, exchange interactions are cal- 
culated for the self-consistent potential and the correct 
Fermi energy, and then also for the same potential but 
assuming the Fermi energy to be slightly higher/lower 
to simulate holes/electrons. The Fermi energy coming 
out from the self-consistent calculation for the equilib- 
rium lattice parameter is simply moved up or down by 
0.136/0.68 eV (0.01/0.05 Ry). In Fig Hi] we show this 
change for both ZB and RS structure Geo.75Mno.25Te and 
A(E F ) = ±0.68 eV (0.05 Ry). In both cases, i.e. Fermi 
energy changed by ±0.68 eV, there is an increase in the 
Mn-Mn interaction. In Fig[12] we compare the results of 
changing he Fermi energy by 0.136 eV and 0.68 eV in 
RS Geo.75Mno.25Te for two different lattice parameters. 
This indicates that both electron- and hole-doping of the 
system would be an efficient way to drive it toward ferro- 
magnetism. Of course, hole- and electron-doping can be 
simulated in other ways as well, for example, by chang- 
ing the effective valence of the atoms. This can be done 
by changing the valence of the magnetic or the nonmag- 
netic atoms. In Fig [13] we show the effect of altering the 
valence of Mn atoms from 7 to 6.8, 6.9 and 7.1 for dif- 
ferent lattice parameters. This is equivalent to studying 
the effect of alloying of Mn with other elements to the 
left or right in the same row of the periodic table using 
virtual crystal approximation (VCA). These calculations 
are self-consistent, as opposed to the above results with 
rigidly shifted Fermi levels, but may suffer from the weak- 
nesses of VCA. We have also included the case where Mn 
concentration in the Ge-sublattice is reduced from 25% 
to 22% and the remaining 3% of the sites on this sub- 
lattice are left vacant. In all these calculations it is the 
magnetic atom Mn, whose concentration is effectively al- 
tered to create the electrons or holes. Next, we explore 
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FIG. 11: (Color online) Change (non self-consistent 
result) in Mn-Mn exchange interactions in RS and ZB 
Geo.75Mno.25Te computed by moving the Fermi energy 

(rigid shift) by 0.68 eV (0.05 Ry) up or down with 
respect to the result of the self-consistent calculation for 
the alloy. Solid lines with solid circles are the results of 
self-consistent calculations for the eauilibrium lattice 
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FIG. 12: (Color online) Change (non self-consistent 
result) in Mn-Mn exchange interactions in RS 
Geo.75Mno.25Te computed by moving the Fermi energy 
(rigid shift) by 0.136 eV (0.01 Ry) and 0.68 eV (0.05 

Ry) up or down with respect to the result of the 
self-consistent calculation for the alloy. Results are 

shown for two different lattice parameters: (a) 
equilibrium lattice parameter 5.964 A, (b) expanded 
lattice parameter 6.125 A. 



the effect of creating holes or electrons by keeping the 
concentration of Mn atoms fixed at 25%. 

We study separately the effect of introducing holes by 
doping the Ge-sublattice and the Te-sublattice. To cre- 
ate holes in the Ge-sublattice, we replace some of the Ge- 
atoms with Cu, and to create holes in the Te-sublattice 
we replace some Te-atoms with Sn or simply vacancies. 
Results are shown in Figs. Q3] and fTS] Fig. [TH shows 



Distance (a) 

FIG. 13: (Color online) Change in Mn-Mn exchange 
interactions brought about by change in carrier 
concentration in various different ways in RS 
Geo.75Mno.25Te. (a) Self-consistent result for 
equilibrium lattice parameter 5.964A, (b) same as (a), 
but by assigning 6.9 electrons to Mn atoms, thus 
simulating hole-doping, (c) same as (b), but by 
assigning 6.8 electrons to Mn atoms, thus simulating 
larger hole concentration, (d) same as (b) or (c), except 

that Mn valence has been moved to 7.1, simulating 
electron-doping, (e) self-consistent result for equilibrium 
lattice parameter, but by replacing 3% of the atoms in 
the Ge-Mn sublattice by empty spheres (vacancies). 



the effect of replacing some of the Ge atoms with Cu, 
thus creating holes. Results show the effect of increas- 
ing hole-concentration as well as the effect of changing 
the lattice parameter or the volume per atom. We have 
considered only the RS all the available experi- 

mental results are for this structure. It is clear that even 
a reasonably low level of Cu-doping of the Ge-sublattice 
can substantially increase the Mn-Mn exchange interac- 
tion, changing it from negative to positive and explaining 
the observed ferromagnetism of Gei-^Mn^Te thin films. 
Fig. [TBI shows the effect of creating holes via doping the 
Te-sublattice with Sn or creating vacancies in this sublat- 
tice. Again, there is a change in the Mn-Mn interaction 
in the direction of ferromagnetism (decreasing antiferro- 
magnetic interactions) . 



VI. SUMMARY OF RESULTS AND 
CONCLUSIONS 

We have examined theoretically the possibility of half- 
metallic ferromagnetism in bulk Gei-^TM^Te alloys, 
with TM being the transition metals V, Cr and Mn. 
FP-LAPW calculations reveal the possibility of half- 
metallicity for some of the ordered alloys. For these 
cases we have calculated the exchange interactions and 
the Curie temperatures. The later calculations apply to 
cases where the Ge and the TM atoms occupy randomly 
the Ge-sublattice. The effect of this disorder is taken 
into account using the CPA. Ferromagnetic interactions 
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FIG. 14: (Color online) Mn-Mn exchange interactions 
in RS Ge-Mn- Te, obtained by putting 2% and 5% Cu in 
the Ge-Mn sublattice. Cu-atoms randomly replace the 
Ge-atoms, thus keeping the Mn concentration fixed at 
25%. The legends (a)-(d) refer to different lattice 
parameters: (a) equilibrium Geo.75Mno.25Te lattice 
parameter 5.964 A, (b) expanded lattice parameter 
6.125 A, (c) and (d) contracted lattice parameters 5.750 
A and 5.534 A, respectively. 



are strongest in the Cr-doped ZB GeTe. The Cr-doping 
case is most promising, as in both RS and ZB struc- 
tures FM interactions dominate, and these remain FM 
with changes in the lattice parameter. The next promis- 
ing case is V-doping in the ZB structure, the RS coun- 
terpart showing strong antiferromagnetic (AFM) inter- 
actions, particularly at and around the equilibrium lat- 
tice parameter. These AFM interactions weaken on both 
sides of the equilibrium lattice parameter value and be- 
come FM only at much higher volume. Our calculations 
for TM=Mn and £=0.25 shows the substance to be AFM 
(for the ordered as well as disordered compounds, with 
the gap being narrower in the disordered case) . We show 
that this AFM behavior is linked to the substance being 
a narrow/zero gap semiconductor/semimetal at this Mn 
concentration. We further establish that the presence 
of uncompensated carriers should drive the material to- 
ward ferromagnetism. The origin of these carriers could 
be vacancies, impurity atoms, and/or structural imper- 
fections. In fact both Ge-Te and Gei-^TM^Te samples 
are known to have large number of holes (p-type carriers) 
due to cation(Ge) vacancies j 15 i 44 ~ — A comparison of the 
results in Figs. H4l and [TBI clearly shows that holes in the 
Ge-Mn sublattice is much more effective in driving the 
system to ferromagnetism than the holes in Te-sublattice. 
This seems to be consistent with the experimental obser- 
vation that all ferromagnetic films of Gei-^TM^Te have 
a large number of temperature-independent p-type carri- 
ers, associated with vacancies in the Ge sublattice. The 
calculated exchange interactions due to carriers originat- 
ing from these causes are sufficient to account for the 
observed ferromagnetism in thin films of Gei-^Mn^Te. 
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FIG. 15: (Color online) Mn-Mn exchange interactions 

in RS Geo^Mno^Te-EZi-a;, where Z stands for 
Sn-atoms or vacancies. Sn-atoms or vacancies occupy 
random positions solely in the Te-sublattice. Legends: 
(a) a;=1.0, a= 5.964 A, (b) £=1.0, a=6.286 A, (c) 
x=0.95, a= 5.964 A, (d) x=0.95, a=6.286 A. 



Our results for the exchange interactions are consistent 
with experiments in that the highest T c reported so far 
seem to be for thin films of Gei-^Cr^Te. Some differ- 
ences between our calculated results and the experiments 
could be ascribed to the fact that our calculations are for 
bulk alloys, while experimental results are mostly for thin 
films. For example, distortions with respect to the bulk 
structures may be present in thin films, in addition to 
the fact that the lattice parameters of the thin films are 
usually larger than their bulk counterparts.— 
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